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Though it is widely taken for granted that it makes sense to separately discuss quark and gluon 
jets, normal jet algorithms lead to a net parton-level jet flavour that is infrared (IR) unsafe. 
This writeup illustrates the problem and explains how the kt algorithm can be modified to 
provide an IR safe parton-level fiavour. As well as being of use in theoretical calculations that 
require a projection of higher-order contributions onto a flavour-channel of the lowest order, 
jet-flavour algorithms also open up the prospect of significant improvements in the accuracy 
of heavy-quark jet predictions. 



1 Introduction 

Over 350 articles on SPIRES refer in their title to "quark-jet(s)" or "gluon-jet(s)" . This pre- 
supposes that such a distinction can be made sensible at parton level. 

It is well known that there is no unique way of defining jets — for example the mapping of 
an n -|- 1 parton configuration onto n-jets (and the decision of whether that mapping should be 
carried out in the first place) becomes ambiguous when all re -|- 1 particles are hard and widely 
separated in angle. This ambiguity persists when trying to identify re flavoured jets from n -|- 1 
partons. Nevertheless one might hope that, if one identifies the flavour of a jet as the sum of 
flavours of its constituent partons, then that flavour will be meaningful, i.e. infrared (IR) safe, 
just as the energies and angles of the jets are IR safe. 

When mapping n-|-l partons onto re jets, the IR safety of the flavour holds trivially in normal 
jet algorithms. However with re -|- 2 or more partons, there can be configurations (fig. ^ with an 
extra large-angle soft qq pair stemming from the branching of a soft gluon, such that the quark 
is clustered into one jet while the anti-quark is clustered into another. Both those jets have 
their flavours 'contaminated' by the soft (anti)quark. Because there is a soft divergence associ- 
ated with the production of the gluon that branched to the large-angle qq pair, a perturbative 



Figure 1: 5 parton configuration clustered to 3 jets, where a large-angle soft qq pair (kijk^) contaminates the 

flavour of two of the jets. 



calculation of thejet flavours will lead to a result that diverges in the infrared. This happens 
notably for cone,^ and Cambridgel^ type jet algorithms. 

The jet algorithm that offers the most scope for resolving this problem is the kt algorithm. It 
introduces a distance measure yij between each pair of particles i and j, recombines the pair of 
particles that are closest, and then repeats the procedure until one has a predetermined number 
of jets (or all jets are sufficiently widely separated). The kt distance measure (in its form for 
e+e~ collisions), is 



2min(K2, 

y"^ = ^^^(l-cos%), (1) 



where Ei is the energy of particle i, 6ij is the angle between particles i and j and Q is the 
centre of mass energy. This choice of distance measure can be justified on the grounds that the 
emission of a gluon is associated with two divergences: one as the gluon energy vanishes, and 
another as the gluon becomes close in angle to any other particle, for example: 
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Such a structure involving two divergences holds only for gluon production. Quark production 
in contrast has only a collinear divergence: 



(note the "max" in the denominator). When one is interested mainly in the kinematics of jets 
this is largely irrelevant, since most of the branchings in an event produce gluons. However when 
investigating flavour, a problem arises because the y^* makes it easy for a soft quark to recombine 
with a hard particle even though there is no corresponding divergence for the production of that 
soft quark. 

A simple solution to the problem is to modify the jet distance measure to reflect the structure 
of divergences in the theory. So we introduce a new "flavour distance" , 

p 2{1 — cos 9ij) J max(£'?, £'?) , softer of i,j is flavoured, , , 

~ Q2 1 mm{Ef,Ej), softer of i,j is flavourless. ^ ^ 

Note that this requires information on the flavour of each object, and so can only be used at 
parton level (an exception is for heavy flavour to which we return later). Figure 121 illustrates 
the Uij distances between various partons in an event with the two distance measures and in 
particular shows how a soft quark has a large yfj with all hard particles in the event. Thus it 
will first recombine with the antiquark of similar softness. This will produce a gluon-flavoured 
object which can recombine with hard objects without changing their flavour. So a modification 
of the distance measure to better reflect the divergences in the theory allows one to obtain an 
infrared-safe definition of jet flavour^ 




Figure 2: Representation of yij distances between various particles in an e'^e event with a kt distance measure 
and the flavour distance measure. Thick lines indicate a large ytj, thin lines a small yij. 



The IR safety can be concretely illustrated by taking a fixed-order parton-level e~^e~ event, 
clustering it to two jets and examining the cross section for cases where the flavour of the two 
jets is not simply that of q, q. If this cross section is plotted as a function of y^^^ the ycut 
resolution threshold in the kt algorithm above which the event is clustered to two jets, the cross 
section should vanish as 1/3* — 0, i.e. as the soft/collinear limit is approached. At order a^, 
fig. this is seen to happen for the flavour algorithms (which actually form a class defined by 
max{Ef,E]) [max{Ei, Ej)]" ■ [mm{Ei, Ej)]^''^ , < a < 2), but not the h algorithm. 
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Figure 3: Proportions of events whose flavour is misidentifled after clustering with various algorithms, as a function 
of the softness of the event: (a) in e^e~, coefficient of O (o^s), calculated with EVENT2'^ (b) in qq qq LHC 

events (with ~ 1 TeV jets) calculated using HerwiJP 

So far we have examined the flavour algorithm just for e^e~ collisions. In hadron-hadron 
collisions (and some DIS contexts) a longitudinally invariant algorithm is used, in which the 
distance measure d^- is obtained by removing the normalisation, taking Ei —>■ ku and 2(1 — 



ti- 



COS 9ij) ^Vij + ^Hj-: t)y introducing an additional beam distance measure diB = k- 
For the flavour algorithm, the kinematic modifications are identical, however the beam distance 
measure that is introduced is more complex: 



(F) ^ I max(/c2.,/c2^(77,)) , 
1 ^Hkti,k^BiVi)) , 



i is flavoured, 
i is flavourless 



(5) 



where one uses a beam hardness as a function of rapidity, k^^{7]), defined in detail in reference^. 
It is more complex to illustrate the IR safety in hadron-hadron collisions than in e~^e~ because 
hadron-collider fixed-order (NLO) programs do not currently provide any information on parton 
flavour. So instead we take HERWIG parton shower events, and look (fig. Ob) at the proportion 
of events in which the reconstructed jet flavours fail to correspond to the original 2^2 event 



jet flavours before showering, again as a function of the event softness. Ah algorithms show 
a failure rate that vanishes as 2/3 —> because in that limit all parton-showering is forbidden, 
however the much faster vanishing for the flavour algorithms is a sign of their IR safety. One 
can also impose hlandness of recombination^ i.e. forbid multi-flavoured recombinations {e.g. 
uu, ud). This has no effect on the 2/3 dependence of the failure rate, but improves the overall 
normalisation. 

There are various applications of jet flavour algorithms. A number of theory calculations 
{e.g. CKKW matchingl^l or hadron-collider resummations^^ have a need to project some high 
fixed-order calculation onto a flavour channel of lower order so as to match with an appropriate 
parton-shower or resummation based on that lower order. 

At hadron level the use of the flavour algorithm is hampered by the need to know the flavour 
of every object, making it inappropriate for experimental discrimination of quark-like versus 
gluon-like hadron-level jets. However there does exist one important hadron-level application: 
heavy-flavour jets. It is feasible experimentally^^to identify all heavy-flavour hadrons in an event. 
Therefore one can use the flavour algorithm treating only heavy-flavour objects as flavoured. The 
infrared safety of the algorithm means that the heavy-flavoured jet distribution will then be free 
of any logarithms of Pt / rrin (other than those in the heavy- flavour parton distribution) , where Pt 
is the jet transverse momentum and uih the heavy-hadron mass. Furthermore, neglecting power 
suppressed terms ~ rrifj/Pf, it can be calculated using a normal light-flavour NLO program*!^ 
(extended to provide access to the flavour information) rather than a dedicated heavy-flavour 
program^^ 

The fact that the new algorithms for safe heavy- flavour jets lead to far fewer logs of Pt/rriH 
than currently used heavy-jet definitions^^ and that light-flavour NLO programs allow the 
remaining logs to be resummed in the PDFs, suggests that heavy-flavour jet distributions can be 
predicted significantly better for the new algorithms, probably reducing the 40 — 60% uncertainty 
typical of current heavy flavour NLO calculations down to the 10 — 20% typical of inclusive light- 
jet distributions at NLO. 
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